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A numerical study of Goldstone-mode effects and scaling functions of the 
three-dimensional 0(2) model * 

J. Engels a , S. Holtmann a , T. Mendes a , and T. Schulzc a 

a Fakultat fur Physik, Universitat Bielefeld, 
D-33615 Bielefeld, Germany 

We investigate numerically the three-dimensional 0(2) model on 8 3 — 160 3 lattices as a function of the magnetic 
field H . In the low-temperature phase we verify the //-dependence of the magnetization M induced by the 
Goldstone modes and determine M in the thermodynamic limit on the coexistence line both by extrapolation 
and by chiral perturbation theory. We compute two critical amplitudes from the scaling behaviours on the 
coexistence line and on the critical line. In both cases we find negative corrections to scaling. With additional 
high temperature data we calculate the scaling function and show that it has a smaller slope than that of the 
0(4) model. For future tests of QCD lattice data we study as well finite-size-scaling functions. 



1. INTRODUCTION 

O(N) models are of general relevance to con- 
densed matter physics and to quantum field the- 
ory, because many physical systems exhibit a 
second-order phase transition with the same uni- 
versal properties. Due to the existence of mass- 
less Goldstone modes in O(N) models with N > 1 
and dimension d = 3 and 4 singularities are ex- 
pected on the whole coexistence line T < T c , H = 
0, in addition to the known critical behaviour at 
T c . Recently these predictions have been con- 
firmed by simulations of the 3d 0(4) model [^). 
A further motivation for studying 3d O(N) mod- 
els is their relation to quantum chromodynamics 
(QCD). The QCD chiral phase transition for two 
light-quark flavors is supposed to be of second or- 
der in the continuum limit and to be in the same 
universality class as the 3d 0(4) model ||. In 
the staggered formulation of QCD on the lattice 
a part of chiral symmetry is remaining and that 
is 0(2). For the comparison to QCD lattice data 
it is therefore important to know the 0(2) and 
0(4) universal scaling functions. 

The 0(2)-invariant nonlinear c-model (or XY 
model) on a d— dimensional hypercubic lattice is 



defined by 

/3W = -J S *- S i - H 



(1) 



<i,3> 



Si is an 2-component unit vector at site i with 
a longitudinal (parallel to the magnetic field H) 
and a transverse component 

S, = Sfu + Si . (2) 

The order parameter of the system, the magneti- 
zation M, is given by 
1 



w =< v^ s 
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= <5« 
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(3) 



There is a longitudinal and a transverse suscepti- 
bility 



XL = H = V« S" 2 > — M 2 ) , 



(4) 
(5) 
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XT = V<S^ 2 >=§. 

In the broken phase (T < T c ) the magnetization 
attains a finite value M (T, 0) at H — 0. Conse- 
quently the transverse susceptibility diverges as 
H- 1 when H -> for all T < T c . It is non-trivial 
that also the longitudinal susceptibility is diverg- 
ing on the coexistence line for 2 < d < 4. The 
predicted Q divergence for d = 3 is 

X l(T<T c ,H) ~ if" 1 / 2 , (6) 
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which is equivalent to an i? 1//2 -behaviour of the 
magnetization near the coexistence curve 

M(T < T C ,H) = M(T,0) + cH 1/2 . (7) 

In finite volumes and H — > the Goldstone modes 
induce strong finite-size effects at all T < T c . 

2. NUMERICAL RESULTS 

Our simulations were done on three-dimensio- 
nal lattices with periodic boundary conditions 
and linear extensions L — 8 — 160 using the clus- 
ter algorithm. In order to eliminate finite-size 
effects we simulated for increasingly larger values 
of L at fixed values of J = 1/T (i.e. fixed tem- 
perature T) and H. For 1/T C we took the value 
J c = 0.454165(4) from Ref. §. 

In Fig. [I] we show the data for the magneti- 
zation as a function of H 1 / 2 for six fixed values 
of J in the low-temperature phase. The picture 
is rather similar to the one obtained in 0(4) 0: 
strong finite-size effects appear for small H and 
persist as one moves away from T c , the results 
from the largest lattices are at first sight linear in 
-ff 1 / 2 , as predicted by Eq. (?]. Very close to H = 
the fixed temperature curves become slightly flat- 
ter, leading to a higher value for M(T, 0) than 
expected from the data at larger H values. This 
behaviour is more pronounced close to T c than at 
lower temperatures. In order to extrapolate the 
data to H — > and V — * oo we apply two dif- 
ferent strategies. The first is to extend the linear 
form in H 1 / 2 , Eq. [?], to a quadratic one 

M(T < T c , H) = M(T, 0) + ci H 1 ' 2 + c 2 H, (8) 

and to fit the data from the largest lattices, which 
we assume to represent data on an infinite vol- 
ume lattice, to this form. The second way to find 
M (T, 0) is just opposite to the first. Here we ex- 
ploit the L or volume dependence at fixed J and 
fixed small H to determine via chiral perturba- 
tion theory (CPT) || the magnetization £ of the 
continuum theory for V — > oo. H = 0, which is 
related to M(T, 0) by 

M(T,0) = A (9 ) 
We observe in Fig. Ill a remarkable coincidence of 
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Figure 1. The magnetization vs. H 1 / 2 in the low- 
temperature region for fixed J and various L. 



the fits according to Eq. || (dashed lines) with 
the CPT results at H = (filled circles). In 
the neighbourhood of T c the results for M(T, 0) 
should show the usual critical behaviour. Since 
we expect here sizeable corrections to scaling H 
we make the following ansatz to determine the 
critical amplitude B (i =T C — T) 

M(T<T c ,0) = Bi0[l + bii" v + b 2 t\ . (10) 

Here and in the following we use the critical ex- 
ponents from Ref. 
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Figure 2. The magnetization at H = vs. T — T c 
with the fit ([h]) (solid line) and its leading part 
(dashed line). 



(3 = 0.3490(6), v = 0.6723(11), u = 0.79(2) .(11) 

A fit to all extrapolated points gives 

B = 0.945(5), h = -0.053(23), b 2 = -0.098(23). 

In Fig. |2| we show this fit and also the leading 
term separately. As the critical point is reached 
the -ff-dependence of the magnetization changes 
to satisfy critical scaling. We therefore fit the 
data from the largest lattice sizes at T c to the 
form 



M(T C ,H) = d c H l/s [l + d l c H^-\ 



(12) 



A further term proportional to H is unnecessary 
here, because the corrections to scaling are much 
smaller than on the coexistence line. The largest 
L data can be fitted very well with the ansatz 
( fi2| ) and the critical exponents as input. We find 
d c = 0.978(2) and d\ = -0.075(5), that is again 
negative corrections to scaling. 

3. THE SCALING FUNCTION 

In the thermodynamic limit, the dependence of 
the magnetization on temperature and magnetic 
field can be expressed || in the form 



Figure 3. The scaling function fc (solid line) and 
results for M/h 1 / 8 at J = 0.55,0.50,0.47,0.465 
and 0.462 (dashed lines), starting with the lowest 
curve. The circles are single data points, the inset 
shows f G for 0(2) and 0(4). 



M = h 1/5 f G {t/h 1/l3S ) 



(13) 



where fa is a universal scaling function and t and 
h are the normalized reduced temperature t = 
(T - T c )/T„ and magnetic field h = H/H . Here 

H = d- s = 1.11(1) , T = B-W = 1.18(2) ,(14) 

which implies /g(0) = 1 and fo{t < 0, h — ► 0) — > 
(— tY^h^ 1 / 5 . Obviously, fa does not account for 
possible corrections to scaling and is the leading 
term in the Taylor expansion in h UUc of a more 
general form 

M/r 1 / 5 = *(thr x l fiS ,h u ' v ') . (15) 

We therefore perform quadratic fits to our data 
in /i wl/c at fixed values of thr 1 /^ 8 in the low- 
tcmpcrature region, where the corrections are 
strong. In the high-temperature region the data 
scale directly. In Fig. |^ we have plotted Mhr 1 ! 8 
from data with H < 0.0075 and 0.43 < J < 0.55 
and the final result for the scaling function Jq. 
An alternative scaling form is that of Widom and 
Griffiths. It is discussed in Ref. |§], where more 
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M = L- f} ^Q 0z (hL 1 /^) + 



(17) 
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Figure 4. (a) M L$l v from reweighted data from 
lattices with different L (solid lines) and the scal- 
ing function Qo (dashed line) at z = 0. (b) Com- 
parison of In Qo to its asymptotic value (line) . 



details of our calculations can be found. 

In Ref. H staggered lattice QCD data for N T — 
4 were compared to the 0(4) scaling function. 
The test failed because the data were indicating 
a steeper scaling function. Since the 0(2) scaling 
function is even flatter than the one for 0(4), as 
can be seen from the inset of Fig. ||, the situa- 
tion will be worse there. A way out may be the 
comparison to finite-size-scaling functions, since 
lattice QCD is presumably far from the thermo- 
dynamic limit. 

3.1. Finite-size-scaling functions 

The general form of the finite-size-scaling func- 
tion for the magnetization is given by 

L-W v $(tL l > v ,hL l /"%L- w ) . 



M 



(16) 

On lines of fixed z = th^ 1 ^ 5 and after expanding 
in we have 



where Qq z is a universal function. Examples of 
lines of fixed z arc the critical line where z = 
and the pseudocritical line, the line of maximum 
positions of the susceptibility xl in the (t, h)- 
planc for V — ► oo. The scaling function contains 
also all information about xl, because 

xl = — = —mr I /cW - -s!g( z ) ■ ( 18 ) 



H 6 



fixed h and 
, of the func- 



Evidently, the maximum of xl at 
varying t is at the maximum point z t 
tion in the brackets of Eq. (|l^). z p is again a 
universal quantity and we find z p = 1.556 ± 0.10 
from (18). As a check we have also determined 
the peak positions on lattices with L = 24 — 96 
which extrapolate to z p — 1.65±0.10 for L — > oo. 
In Fig. ||a we show ML 13 / 11 from the reweightcd 
data on the critical line for various L values and 
also the scaling function Qo(z = 0). For L — > oo 
the scaling function Qq z is related to fa by 

Qo, - faiz^hL 1 /^) 1 ' 6 . (19) 

From Fig. ||b we see that Qo(z = 0) seems to be 
asymptotic in the whole variable range consid- 
ered. At present we investigate the 
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